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AN EXPERIMENTAL INVESTIGATION OF TRANSONIC FLOW PAST 
TWO-DIMENSIONAL WEDGE AND CIRCULAR-ARC SECTIONS 
USING A MACH-ZEHNDER INTERFEROMETER 1 

Bv Arthur Earl Bryson, Jr. 


SUMMARY 

Interferometer measurements are given of the flow fields near 
two-dimensional wedge and circular-arc sections at zero angle 
of attack at high-subsonic and low-supersonic velocities. Both 
subsonic flow with local supersonic zone and supersonic flow 
with detached shock wave have been investigated. Pressure 
distributions and drag coefficients as functions of Mach num ber 
have been obtained. The wedge data are compared with the 
theoretical work on flow past wedge sections of Guderley and 
Yoshihara, Vincenti and Wagoner , and Cole. 

It is shown that the local Mach number at any point on the 
surface of a finite three-dimensional body or an unswept two- 
dimensional body , moving through an infinite fluid , has a 
stationary value at Mach number 1 and , in fact, remains nearly 
constant for a range of speeds below and above Mach number 1. 
On the basis of this concept and the experimental data, pressure 
distributions and drag coefficients for the wedge and circular-arc 
sections are presented throughout the entire transonic range of 
velocities. 

INTRODUCTION 

DIFFICULTIES OF THEORY AND EXPERIMENT IN THE TRANSONIC RANGE 

OF VELOCITIES 

The difficulties inherent in studying transonic How are 
well-known. Theoretical analysis is made difficult by the 
nonlinearity of the differential equations of compressible 
fluid motion. This nonlinearity leads to a change-over in 
type of the differential equations from elliptic to hyper- 
bolic when transition is made from subsonic to supersonic 
speeds. Since the essential feature of transonic flow is this 
mixed subsonic-supersonic character, it is obvious that no 
linearization of the differential equations (at least in the 
physical plane) can adequately describe the flow. 

Wind-tunnel studies in the transonic range are made 
difficult by the large lateral extent of the perturbation flow 
field around bodies in this range. This means that models 
which are small compared with the test section must be 
used. Even then there is still a range of speeds from just 
below M m = 1 to just above M ao = 1 where the model and/or 
its support configuration are “choked,” that is, where local 
supersonic zones embedded in the subsonic field extend from 
the model to the tunnel walls, or, in the supersonic case, 


1 Supersedes NACA TN 2560, “An Experimental Investigation of Transonic Flow past 
Arthur Earl Bryson, Jr., 1951. 


where embedded subsonic zones extend to the tunnel walls, 
or shock waves, reflected from the walls, impinge on the 
model. Some progress has been made recently in modifying 
wind-tunnel test sections so as to minimize these effects, 
but, on the whole, the majority of good test data in the 
range very close to M a = 1 has so far come from free-flight 
tests. Some good transonic data are available, however, 
from transonic-bump tests made in wind tunnels (reference 
1). Using small models usually results in low Reynolds 
numbers so that difficulty is often experienced in extrapo- 
lating data to full-size Reynolds numbers; this seems to be 
particularly true of the transonic speed range since the 
effects of boundary-layer and shock-wave interactions seem 
to be quite large there (references 2 and 3). 

In this paper it is shown that in many instances tests need 
not be made in the region very close to M a> = 1 since the 
flow in this range can be inferred from testing below and 
above this range and using an interpolation based on the 
fact that the local Mach number at any point on the surface 
of unswept two-dimensional bodies and finite three-dimen- 
sional bodies has a stationary value at M O0 = 1. 

EXISTENCE OF POTENTIAL TRANSONIC FLOWS 

Guderley (reference 4) has made a detailed investigation 
of the possibility of smooth transonic flows (i. e., subsonic 
flows with an embedded supersonic zone in which no shock 
waves appear). He proposes that such smooth flows are 
exceptional, that they are discrete cases occurring for only 
particular body shapes at particular free-stream Mach num- 
bers. Any perturbation of the shape with the Mach number 
held constant (or vice versa), Guderley claims, would result 
in a shock appearing in the flow. This bears an analogy to 
the well-known Busemann supersonic biplane which theoret- 
ically has no shocks (and hence no drag) at a discrete value 
of free-stream Mach number and angle of attack (reference 
5, p. 154). Guderley’s proposal is still controversial (e. g., 
see the paper of Sears who has made a critical survey of the 
work to date on the existence of transonic potential flows 
(reference 6)). 

It is obvious that the potential flow must break down for 
a given body shape at some Mach number less than 1. The 
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argument whether this breakdown occurs precisely when a 
supersonic region first appears on the body or at a slightly 
higher Mach number seems somewhat academic (although 
very interesting), since it is well-known experimentally that 
the drag-rise Mach number (i. e., the Mach number where 
noticeable shocks first appear) is very close to the critical 
Mach number (i. e., the Mach number at which sonic velocity 
first appears on the body) for most bodies without surface 
slope discontinuities. 

Kuo (reference 7) proposes that supersonic compression 
is unstable to disturbances; that is, a supersonic region on 
a body in subsonic flow must end in a shock with no com- 
pression occurring in the supersonic flow ahead of the shock. 
There seems to be ample experimental evidence to show 
that this is not strictly true since, for example, the com- 
pression region of a X-shock is clearly supersonic. However, 
the X-shock configuration is believed to be a phenomenon 
associated with laminar-boundary-layer and shock-wave in- 
teraction; with turbulent boundary layer (a condition more 
closely approaching nonviscous flow) hardly any noticeable 
supersonic compression occurs before the shock ending the 
supersonic zone (see reference 8). 

CHOICE OF MODELS 

Two-dimensional flow is much simpler to handle than 
axially symmetric flow both in theoretical work and in inter- 
ferometry. Hence it was decided to study two-dimensional 
flows despite the well-known difficulties in approximating 
two-dimensional flow in a wind tunnel. 

Because of the considerations mentioned previously it 
was decided to test very small models which would be of 
such a shape that viscous influences would not materially 
affect the flow over them. This led to the choice of “half 
airfoils” — wedges and circular-arc sections followed by 
straight sections. These models have favorable pressure 
gradients on their surfaces over most of the transonic 
range so that boundary-layer separation, if it does occur, 
will only occur because of shock-wave influence. Further- 
more, such separation will occur downstream of the part of 
the body being studied and hence will not affect the measure- 
ments. Certain viscous effects will still be evident, however, 
for instance, the effective rounding off of the shoulders and 
leading edges of the wedge models. 

Both theoretical advantages and practical need make the 
study of thin sections desirable. Consequently, the semi- 
wedge angles chosen were 4j°, 7|°, and 10° (a 26.6° wedge 
was also used in order to make a comparison with some avail- 
able theoretical work on a wedge of this angle) . The circular- 
arc section chosen was essentially the front half of an 8.8- 
percent-thick biconvex circular-arc airfoil, followed by a 
straight section. Models of sections much thinner than this, 
with the same chord lengths used, run into structural 
difficulties and also the ratio of boundary-layer thickness to 
model thickness becomes large enough to cause considerable 
deviation from nonviscous flow. 

TRANSONIC-FLOW THEORY AND EXPERIMENTS 

The investigations of Von Karman, Busemann, Guderley, 
Frankl, and many others have contributed significantly to 
methods of approach which can be used to study transonic 


flow (references 9 to 14). The detailed numerical calculations 
for specific cases made by Maccoll and Codd, Emmons, 
Drougge, Drebinger, Guderley and Yoshihara, and Vincenti 
and Wagoner (references 15 to 21) have helped to dispel the 
idea of a “sonic barrier.” Recently Cole at GALCIT has 
given an analysis of the flow past wedge sections at high- 
subsonic speeds (reference 22). By combining the results of 
Guderley and Yoshihara’s, Vincenti and Wagoner’s, and 
Cole’s calculations, the flow past thin wedge sections can 
be given completely through the transonic range permitting 
a comparison with the present experiments. Some of the 
investigations mentioned above will be discussed in more 
detail further on in the present paper. 

Available experiments in the transonic range on thin wedge 
sections are surprisingly few. Pack (reference 23) describes 
some interferometric experiments on 10° and 20° semiangle 
wedges made at Braunschweig. His subsonic data appear 
to be good, but the flow in the supersonic interferograms 
appears to be somewhat nonuniform and not very closely 
two-dimensional; only one supersonic Mach number was 
tested where detached shocks occurred. His conclusion that 
the plp m distributions on the surface of the 20° semiangle 
wedge are very much the same for M 00 = 0.803 and M m = 
1.40 is interesting, but the statement that this agrees with 
the theoretical predictions of Maccoll and Codd is incorrect 
since they indicated that the p/p 0 distributions would be 
nearly the same. 

Griffith at Princeton has just recently published the 
results of some very carefully done experiments on flow past 
wedge sections of semiangles of 7°, 10°, 20°, 30°, 45°, and 
90° (and several other shapes) with detached shock waves 
(reference 24). These experiments were done in a shock 
tube and interferograms are presented of the flow fields. 
The experiments clearly show that the shape of the detached 
shock and its detachment distance from the sonic point on 
a wedge depend only on the body thickness and the Mach 
number (not the wedge angle) when the Mach number is 
well below the shock-attachment Mach number. This is in 
general agreement with Busemann ’s considerations in his 
paper on detached shock waves (reference 10). 

Liepmann, Ashkenas, and Cole (reference 8) made some 
careful pressure measurements on the surfaces of 6- and 
12-percent-thick biconvex circular-arc airfoils at zero angle 
of attack at high-subsonic speeds in connection with studies 
of shock-wave and boundary-layer interaction. Some of the 
results of their tests are combined here with corresponding 
low-supersonic test results from the present investigation to 
indicate the behavior of the pressure distribution on circular- 
arc airfoils at zero angle of attack through the entire transonic 
range. 
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ing out the experimental work. Some of the results of these 
experiments have already been reported in reference 25. 

SYMBOLS 


a 

c 

Cd 

C D 

C v 

c 

1 p 

k 

l 

M 

n 

V 
<1 

t/c 

u 

V 

*, V 


y 

a 

7 

0 

X 


sound velocity 
airfoil chord 
pressure-drag coefficient 


reduced drag coefficient ^^ 5 73 ^ ^) 


pressure coefficient 
reduced pressure coefficient 


/ ( 7 +l) 1/3 fiA 

v m 2/ 3 ) 


Gladstone-Dale constant 
model span 
Mach number 
index of refraction 
pressure 

dynamic pressure 
airfoil thickness ratio 

horizontal component perturbation velocity (per- 
turbation from a*) 

vertical component perturbation velocity 
Cartesian coordinates, origin at leading edge of 
profile 

reduced vertical distance ([( 7 + l)t/c] 1/3 y) 
angle of attack 

ratio of specific heats (1.4 for air) 
semiwedge angle 

wave length of monochromatic light used on inter- 
ferometer 


reduced Mach number 


/ M 2 - 1 

\(y+l) 2,3 (t/cy'\ 


p density 

Subscripts and superscripts: 

( ) ro conditions in free stream 

( ) 0 reservoir conditions 

( )o reservoir conditions behind a shock wave 

( )* conditions at sonic velocity 

Symbols used without subscripts indicate local conditions. 


APPARATUS AND METHODS 

WIND TUNNEL 

The measurements were made in the GALCIT 4- by 
10-inch transonic wind tunnel. For a description of the 
tunnel and the flexible nozzle employed see reference 26. 
The tunnel can be run at both subsonic and low-supersonic 
velocities with continuous Mach number variation through 
use of the flexible nozzle and a variable second-throat nozzle 
downstream of the test section. 

MODELS 

The models used were half airfoils followed by straight 
sections. Four of the models were wedges (semiangles 
4.53°, 7.56°, 10.00°, and 26.57°) followed by straight sections 
and the fifth was half of a biconvex circular-arc airfoil (8.80 
percent thick) followed by a straight section (see fig. 1). 
The distance from the leading edge to the point where the 
straight section began was of the order of % inch for all five 
models. The models were made of tool steel and were very 
carefully machined and lapped so as to give exact cylindrical 


4 . 53 ° .060 

.372 • 


.745 


756 ° -058 



-.219 


.750 



.750 

Figure 1. — Geometry of sections tested. All dimensions are in inches. 


surfaces. Two pressure orifices on opposite sides of the air- 
foil were placed exactly the same distance from the leading 
edge to aid in setting the model to zero angle of attack by 
balancing these pressures on an alcohol U-tube. Because 
of the very short chord lengths vernier-protractor measure- 
ments of the opening angles of the leading edge were of 
doubtful accuracy, so the angles were measured by letting 
the leading edge split a beam of parallel light and measuring 
the position of the reflected spots on a wall behind the 
model. In this manner the angles could be measured to 
±0.03°. 

INTERFEROMETER 

The interferometer used in this investigation is described 
in references 27 and 28. One of the main features of this 
interferometer is that both light beams are passed through 
the test section, one over the model and the other ahead of 
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the model in the uniform flow field, that is, where the velocity 
is nearly the free-stream velocity. The advantages of this 
are: (1) The fringe shifts are in relation to the free-stream 
density and (2) the effects of the side-wall boundary layers 
are approximately canceled out since both beams traverse 
nearly the same boundary la}^er at each side window. This 
leads to improved accuracy when the interferograms are 
evaluated on the basis of the absolute value of the fringe 
shift from no-flow conditions. For these tests finite-fringe 
interferograms were used and another method of evaluation 
was devised which is much simpler and more accurate than 
the above-mentioned technique. Infinite-fringe interfero- 
grams, while they give the constant-density contours im- 
mediately, are less accurate than the superimposed finite- 
fringe interferograms because any optical inaccuracies in the 
system cause the contour fringes to be distorted. These 
inaccuracies are calibrated out in the superimposed finite- 
fringe interferograms. Also there are times when one does 
not know whether the density increment between contours 
of an infinite-fringe interferogram is positive or negative; 
this trouble does not arise with the finite-fringe interfero- 
grams. A typical finite-fringe interferogram is shown in 
figure 2. 



Figure 2.— Typical finite-fringe interferogram. 8.8-percent circular -arc section at =1.200. 

METHOD OF EVALUATION OF INTERFEROGRAMS 

The method of evaluation used here depends on two 
techniques: (a) Photographic superposition of disturbed 
and undisturbed interferograms and (b) fringe identification 
by a pressure measured on the model. 


Direct photographic superposition of a “no-flow” finite- 
fringe interferogram on a “with-fiow” finite-fringe interfero- 
gram gives rise to dashed shadowy lines (the dashes being 
where the dark fringes of one picture cross the light fringes 



Figure 3.— Typical superimposed finite-fringe interferogram. 10° semiangle wedge at 

Af oo =1.278. 


of the other); see figure 3 for an example of this type of 
picture. These shadowy lines can easily be shown to be lines 
of constant density for two-dimensional flow and are the 
same contours as would be obtained on an infinite-fringe 
interferogram made with perfect optical surfaces. The 
increment in density between these shadowy lines is a 
constant dependent only on the span of the model and the 
wave length of the monochromatic light being used. This is 
easily shown since the difference in optical path lengths of 
the light rays between two adjacent constant-density con- 
tours must be 1 wave length of the light being used. For 
two-dimensional flow the difference in optical path length 
will simply be /An , where l is the span of the model and An is 
the difference in index of refraction between the two light 
paths. Thus 

lAn=\ ( 1 ) 

But the relation between index of refraction and density in a 
gas is given by 

n—l=kp (2) 

where k is the Gladstone-Dale constant (a function of the 
i light frequency and type of gas). Therefore 
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A p=\/kl 


( 3 ) 


Ap — \/kl 


where A p is the difference in density between two adjacent 
constant-density contours. For these experiments 


Now in any part of the flow field where the stagnation 
density is constant along a streamline, 


\ = 5461 A (mercury green line) 

&=0.1162 cu ft/slug 
1= 3.50 in. 
so 

Ap/p 0 = 0.0250 per fringe shift 

where p o = 0. 00211 slug per cubic foot was the usual tunnel 
stagnation density. 

The advantage of photographic superposition is not only 
in time saved but also in increased accuracy of evaluation. 
Any slight changes in fringe spacing or fringe orientation 
with respect to the no-flow interferogram which occur before 
the with-flow interferogram is taken can be almost exactly 
canceled out by causing the two superimposed interfero- 
grams to coincide exactly in a region where it is known that 
the flow was uniform, since in such regions there should be 
no isopycnic contours. This is particularly easy to do for 
supersonic flow if a portion of the flow field ahead of the nose 
shock wave is included in the interferogram. For subsonic 
flow care must be taken to include enough of the flow field 
ahead of the model in the interferogram to have some of the 
nearly undisturbed flow field for comparison; this was quite 
simple to do for the small, thin models used in these tests. 

The actual superposition technique used here was first to 
make a print (3)i times enlarged) of the with-flow interfero- 
gram. This print was then placed under the enlarger and 
the no-flow interferogram negative was put into the enlarger. 
By changing the enlargement scale and moving the with-flow 
interferogram under the enlarger the fringes were made to 
coincide exactly in the regions of uniform flow. The con- 
stant-density contours could then be drawn in on the print. 
Alternatively, the first print could be made on transparent 
paper (Ansco Reprolith Ortho was used) and when the super- 
position was accomplished a piece of photosensitive paper 
was slipped under the transparent print and a print of the 
two interferograms was obtained. This was the technique 
used for figure 3. 

In order to identify the density values with the fringes a 
pressure tap was placed on each model approximately half- 
way from the leading edge to the shoulder (a region where 
the pressure gradient was expected to be large). From the 
pressure reading the density at the pressure tap was calcu- 
lated using reservoir fluid properties (taking into account 
entropy changes through shock waves). The pressure tap 
will always lie between two fringe contours or on a contour, 
so that, bj^ knowing the density increment between fringe 
contours, the values of the density on the adjacent contours 
can be obtained by interpolation. The whole interferogram 
is determined once the density is known on one contour (ex- 
cept for the shifts through shock waves). 


Mi+v-r «> 

so 

^=-(l+l=T M 2 ) T_1 MdM (5) 

Hence the increment in Mach number between adjacent 
contour fringes is given approximately by 


7 



since 

A p X 

p 0 klpo 


This function has a minimum at M =^- which is 

M=0.914 for air ( 7 = 1.4). A graph of this function is 
shown in figure 4. ]XTote that the Mach number increment 
per fringe for these tests was always closely equal to 0.05. 



Figure 4.— Increment in Mach number per fringe contour against local Mach number for 
air ( 7 = 1 . 4 ). For these tests X/fc/p„=0.025zfc0.001. 


Similarly, the expression for the increment in pressure 
coefficient between adjacent contour fringes is approxi- 
mately 

x 


AC P 


MJ 1 - 


y— 1 


-A- klp 0 


( 7 ) 


M 2 


r 


For values of M close to M m , this expression has a minimum 




which is 1.832 for air. 


INTERFEROMETER SENSITIVITY 

It is interesting to note that the interferometric method 
has its greatest sensitivity in the transonic range. As 
pointed out previously, the density increment between two 
adjacent contour fringes is a constant 


SIDE-WALL BOUNDARY-LAYER EFFECT ON APPROXIMATING 
TWO-DIMENSIONAL FLOW 

A close approximation to two-dimensional flow over the 
whole span of the model was required since the interferometer 
integrates the value of the density from wall to wall. In a 
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nonviscous fluid letting the model extend from wall to wall 
would theoretically give two-dimensional flow over the whole 
span. If the model did not span the whole tunnel, the flow 
would correspond to that past a model of infinite span with 
periodic gaps in it where the gaps were equal to twice the 
distance from the edge of the model to the wall. The effect 
of the side-wall boundary layers, for a model that does not 
span the tunnel, is roughly to decrease the size of this gap. 
Approximately, the gap size would be decreased by twice 
the displacement thickness of the wall boundary layer. By 
making the gap between the edge of the model and the wall 
approximately equal to the wall-boundary-layer displace- 
ment thickness, one might hope to approximate closely two- 
dimensional flow over the span. This phenomenon is, of 
course, very much more complicated than this, particularly 
in the supersonic case where the shock waves interact with 
the wall boundary layer. However, by taking circular- 
cylinder and wedge models and varying the gap size in in- 
crements of He inch, it was found that the detached bow wave 
became closely two-dimensional when the gap size was / inch 
( i. e., there was no blur ahead of or behind the shock pictures) 
which is almost exactly the bound ary -layer displacement 
thickness when measured without a model in the test section. 
When the gap was He inch the shock was blurred ahead of 
the main shock and when the gap was He inch it was blurred 
behind the main shock. These tests were further sub- 
stantiated by some schlieren pictures, which Mr. Walter G. 
Vincenti of the NACA Ames Aeronautical Laboratory 
kindly made available, showing a view looking down on a 
wedge model so that the leading edge of the detached shock 
appeared as a line; by varying the model span a discrete 
value of the span was found where this line was almost exactly 
parallel to the leading edge of the model, while for just 
slight variations from this gap size the shock was curved 
forward or backward. Figure 2 shows a finite-fringe inter- 
ferogram of the circular-arc section with a detached shock 
where the definition of the shock wave was unusually sharp. 
This is strong, but, of course, not conclusive, evidence that 
the flow was closely two-dimensional over most of the span. 
Further evidence that the flow differed from two-dimensional 
flow only slightly is given in the next section. 

SIDE-WALL BOUNDARY-LAYER EFFECT ON INTERFEROGRAM 
EVALUATIONS 

A result of the method of interferogram evaluation de- 
scribed above is that the effect of the side-wall boundary layer 
is approximately canceled out, since the over-all fringe shift 
from no-flow conditions is unimportant, only the relative 
fringe shifts from a point of known density being used. This 
is strictly true only if the integrated side-wall boundary -layer 
density, defined by 



where y is the direction perpendicular to the tunnel wall and 
y= 0 is the wall, is the same over the entire field of view of 
the interferometer. Obviously, this can never be exactly 
true since the pressure field caused by the model, the 
boundary-layer growth, and the shock-wave and boundary- 
layer interaction all tend to change this value. An indica- 
tion that all these effects might be small was obtained from 


the model tests where pressures were measured at two 
points on the model in the center of the span, where the 
flow is closely two-dimensional; the density increment be- 
tween these two points on the model was compared with 
the density increment given by the interferogram. The 
standard deviation from zero of the difference between these 
two increments over the whole range of test Mach numbers 
was about 1 percent of the stagnation density. Also, the 
values of pressure-drag coefficient obtained interferomet- 
rically for the attached-shock-wave cases checked the 
oblique-shock theory very closely, and it is well-known that 
the oblique-shock theory checks experiment quite well. 

DETERMINATION OF FREE-STREAM MACH NUMBER 

An interesting result of the method of evaluation just 
described is that the free-stream Mach number in subsonic 
flow can be determined from the interferogram and the 
measured pressure on the model, provided a large enough 
field of view ahead of the model is obtained in the inter- 
ferogram. This can be done by noticing that a certain 
number of compression contours appear around the leading 
edge and then expansion contours follow these toward the 
back part of the airfoil; the center fringe corresponding 
to free-stream density can then be traced out into the flow 
field (see, e. g., figs. 9(a) to 9(d) of the 10° wedge in subsonic 
flow). The exact value of the density can be determined on 
this fringe as described previously and, hence, knowing 
the stagnation density in the settling chamber, the effective 
free-stream Mach number can be determined from the 
isentropic-flow relations. It is believed that this effective 
Mach number is a good approximation to the free-flight 
free-stream Mach number and would give the same flow as 
that measured in the wind tunnel for the very small models 
used in these tests. 

This method is more accurate at high-subsonic speeds 
than at low speeds since more contour lines are obtained on 
the airfoil at the higher speeds (see above discussion). The 
estimated accuracy in determining free-stream Mach num- 
ber in this way was ±0.01 for the range of subsonic Mach 
numbers tested. 

The free-stream Mach numbers for the supersonic tests 
were obtained by calibrating the flexible-nozzle jack settings 
against Mach number with a static-pressure probe in the 
center of the tunnel. The probe was traversed upstream 
and downstream in the region where the models were to be 
tested and an average Mach number was obtained there. 
The standard deviations from this average value were of the 
order of ±0.005 in Mach number for the range of supersonic 
Mach numbers tested. 

WIND-TUNNEL CHOKING 

In all the subsonic testing the embedded supersonic zone 
was not allowed to touch the upper or lower walls. In one 
or two of the low-supersonic tests there was a question 
whether the embedded subsonic zone touched the ceiling or 
not. In case it did, it is well-known that in such cases the 
detached shock changes its curvature near the ceiling so as 
to come in nearly normal to the walls. Since the models 
were so very small (He in. thick compared with the 10-in. 
height of the tunnel), it is believed that the effect of this on 
the pressure distribution was negligible. 
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REYNOLDS NUMBER 

The value of the Reynolds number for all of these tests 
was approximately 60,000 based on the chord of the model. 
The boundary layer on the models was laminar and no 
effort was made to trip the boundary layer to make it 
turbulent. The compression region in the shocks shown in 
the high -subsonic-flow int erf programs is believed to be 
associated with the laminar boundary layer, as mentioned 
previously. 

THEORETICAL WORK ON TRANSONIC FLOW 

RELAXATION CALCULATIONS 

In 1946 Maccoll presented a paper at the Sixth Interna- 
tional Congress for Applied Mechanics in which he described 
a relaxation calculation of the compressible flow past a 20° 
semiangle wedge followed by a straight section at Mach 
numbers of 0.7 and 1.5. The flow field in both cases con- 
tained both subsonic and supersonic velocities. His main 
assumptions were: (1) Sonic velocity occurs at the shoulder 
and (2) the streamlines of the flow are perpendicular to the 
sonic line (i. e., the line where sonic velocity occurs in the 
flow). The first assumption can be shown to be correct (see 
reference 22) so that, indeed, it is not an assumption. The 
second assumption, as Maccoll realized, was only approxi- 
mately correct for M m = 1.5 and certainly quite incorrect 
far away from the wedge at M ao = 0.7 (since the assumption 
leads to an infinite supersonic region above the wedge). In 
effect, his solution at M m = 0.7 was “choked” in the sense 
that the back part of the body could have no influence on 
the front part. It is well-known that for bodies at high- 
subsonic speeds a finite, closed supersonic region occurs in 
the flow, so that the sonic line makes all angles possible with 
the streamlines, including 0°. The method of solution used 
was to assume positions of the shock wave and sonic line, 
calculate the residues in the relaxation net using the isen- 
tropic-flow equations (an approximation since flow behind a 
curved shock is not isentropic), and then readjust the shock- 
wave and sonic-line location, calculate again, and so forth, 
iterating until the solution closely repeated itself. Maccoll 
found that the 'plpo distribution on the wedge surface at 
M m = 1.5 was nearly identical with the p/p 0 distribution at 
M m = 0.7. This led him to propose that the pressure in the 
transonic region, on bodies with distinct corners, varied as 
the stagnation pressure and he presented a drag curve 
through M m = 1 for the 20° semiangle wedge calculated on 
this basis. 

Drougge in 1948, following Maccoll, calculated the flow 
past a finite cone of 45° semiangle with detached shock wave 
at M.^1.80 and M O0 = 2A5 ) using the same assumptions 
as Maccoll (reference 18). He also made experiments on 
this cone and found the agreement with his theory rather 
good. He made several tests at lower supersonic Mach 
numbers also and found that the plp 0 ' distribution on the 
cone surface did remain nearly constant except as the Mach 
number became close to the attachment Mach number. 

Drebinger in 1950 showed how to calculate, by relaxation 
techniques, the flow past finite cones and wedges with de- 
tached shocks, eliminating the isentropic-flow assumption 
and the assumption on the streamlines being perpendicular 
to the sonic line (reference 19). He calculated a specific 

243695—53 2 


example — a 26.6° semiangle wedge at M m = 1.440 — and 
checked the calculated shock-wave shape and position exper- 
imentally. His calculations showed that, even for the de- 
tached-shock case, the streamlines differed from being per- 
pendicular to the sonic lines by angles as large as 30°. His 
calculation was checked in detail experimentally in these 
tests and agreement was found to be excellent. 

TRANSONIC PERTURBATION THEORY 

By assuming that the velocity component parallel to the 
free-stream direction differs only by a small quantity u from 
a *, the critical velocity, and keeping only the highest-order 
terms in the differential equation, the equations of two- 
dimensional irrotational fluid motion are reduced to 


f du.bv J) 

(7+1) ^^ + ^ =0 ( 

bu bv . ( 

by bx~~ ) 


(9) 


It was from these equations that Von Karman and Guderley 
independently arrived at the transonic similarity laws 
(references 9 and 4). For two-dimensional steady flow past 
sections whose shape functions are the same, these laws 
imply that 


M 2 -l ( MJ - 1 ) 

[( 7 +l)*/c ] 2/3 J \[(y+l)tm 


( 10 ) 


where M is the local Mach number on the surface of the 
section. The similarity in pressure and drag coefficients 
is then 


(7+l) 1/3 ^ A MJ - 1 ) 
m m J \[(y +i)t/c]*'A 

(7+1 y /z c D f mj- i ) 
m b,z (Kt+1 y/cH 


(11) 

( 12 ) 


These quantities will be called reduced local Mach number, 
reduced free-stream Mach number, reduced pressure coef- 
ficient, and reduced drag coefficient, respectively, using 

symbols £, £„, C v , and C D . 

By interchanging dependent and independent variables 
in the perturbation equations, the problem becomes linear: 


where 


_ bg , d| =0 

bv du 


d£_d]/ = r l 

bv bu 


\ 

i 


m=(y+1) 


u 

T* 


(13) 


fl = (7+l) 


and, by eliminating x by differentiation, the Tricomi equa- 
tion is obtained: 


_ b 2 y b 2 y n 

U ^-= 2 +^-= 2=0 

bv 2 bu 2 


(14) 


The main difficulties with this liodograph ( u , v) plane are: 
(a) The mapping of physical boundaries into the liodograph 
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plane is, in general, not known until the solution to the 
problem is known so that it is not known where to apply the 
boundary conditions in the hodograph plane and (b) the 
mapping is often multivalued, complicating the solution. 
Two interesting cases are known where these difficulties 
are avoided. They are: (a) The free jet, studied by Tscliap- 
lygin in 1905, and (b) the finite wedge, studied recently by 
Guderley and Yoshihara, Vincenti and Wagoner, and Cole. 
These latter studies came to the author’s attention after the 
present experimental study of the finite wedge in transonic 
flow had begun and served to make the study more interest- 
ing since the data could then be compared with the theo- 
retical results. 

THEORETICAL STUDIES OF TRANSONIC FLOW PAST THIN WEDGE SECTIONS 

Guderley was the first to formulate the problem of the 
thin finite wedge in the hodograph; he and Yoshihara found 
an approximate solution to the problem of the flow past a 
thin double-wedge profile at zero angle of attack at Mach 
number 1 using the transonic perturbation equations 
(reference 20). 

Vincenti and Wagoner considered the thin double-wedge 
profile at zero angle of attack for low-supersonic Mach 
numbers where the shock wave is detached (reference 21). 
Their solutions were effected by relaxation calculations in 
the hodograph plane. Here the bow shock wave and the 
sonic line are fixed boundaries (their positions are not 
known originally in the physical plane) and the boundary 
condition on the shock is the slope of the streamlines (or 
the lines y= Constant). This boundary condition was 
first shown by Busemann, who aptly called the configura- 
tion a “hedge hog.” 

Cole (reference 22) has recently given a simple approxi- 
mate analytical solution to the flow past a thin symmetrical 
wedge followed by a straight section at high-subsonic 
speeds (M m ^l). His solution satisfies the Tricomi equa- 
tion and the boundary conditions on the wedge and at infinity 
but not the boundary conditions on the sonic line. Effec- 
tively, his solution gives a finite vertical sonic line from the 
shoulder which is also a limiting line. Cole has indicated 
that this solution is the singular part of the solution in the 
hodograph and as such is most likely the main part of the 
solution. It is interesting to note that the drag-curve slope 
and curvature at M m — 1 obtained from Cole’s solution 
agree exactly with the values obtained from the simple 
physical considerations of the next section. Also, the pres- 
sure distribution on the wedge at M m = 1 agrees within 1 
or 2 percent with that obtained by Guderley and Yoshihara. 

Since the back half of a double-wedge profile has only a 
very weak influence on the pressure distribution on the 
front half for M 00 > 1 (only through the “last Mach wave” 
from the shoulder point to the sonic point on the detached 
shock), it is reasonable to take the solution of the double 
wedge at M m ^ 1 and use the front-half solutions in con- 
nection with Cole’s results for M m ^ 1 for the wedge fol- 
lowed by a straight section and thus have a solution for the 
latter semi-infinite body completely through the transonic 
range. By using linearized subsonic theory and the shock- 


expansion supersonic theory, the zero-angle-of-attack flow 
is obtained for all possible values of M m . 

Tsien and Baron (reference 29) have shown that the 
shock-expansion theory can be expressed in the transonic 
similarity form for thin bodies in pure supersonic flow near 
M od = 1. 

Von K&rmdn (reference 9) has indicated also how linear- 
ized subsonic- and supersonic-flow results may be written 
in the transonic similarity form since, from the Prandtl- 
Glauert similarity, in linearized subsonic theory, 

C„= - t/C 1 — M J~) (15) 

-y l — MJ J \ c ej 

and, in linearized supersonic theory, 



and from the expressions for reduced pressure coefficient 

( 7 + 1 ) 1/3 

and Mach number, multiplying both sides by - yj c y/ z ~ > 
these equations may also be written as: 


(7+1 ) 1I3 C P 
0 t/c ) 2/ 3 


(7+1 Y IZ C P 
(i t/c) 2/z 


but 


[(7+l)/+] 2/3 

V 1 -MJ 


A7^W+wM. h + 1) i Hi (150) 

/ [(7+1)</c] 2/3 
V MJ - 1 


|W[(vhwM (7+1) 3 ?! (16a; 


* (t+i) ,/3 <?„ 

p (*/c ) 2/3 


L 


M n 


1 


-[( 7 +1)*/c] 2/3 

y=[(y+l)t/c] m y 

so equations (15) and (16) may be written in transonic form 


c-=r( H-V=r.) 

(17) 


(18) 


The subsonic pressure-distribution and drag-coefficient 
curves have been calculated here from Cole’s analytical 
expressions and, combined with the results of Guderley and 
Yoshihara, Vincenti and Wagoner, and Tsien and Baron, 
the curves for reduced pressure and Mach number distribu- 
tion and reduced drag coefficient 2 are given in figures 5 to 8 
for the finite wedge followed by a straight section. 

1 The reduced drag coefficient given in figure 8 is that for the half wedge and is equal to 

-+ (t) 
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Figure 5.— Theoretical reduced local Mach number distributions on a wedge near Mach 
number 1. Data from references 20 to 22 and 29. 


Figure 7. — Theoretical reduced pressure-coefficient distributions on a wedge near Mach 
number 1. Data from references 20 to 22 and 29. 




Figure 8.— Theoretical reduced drag coefficient against reduced Mach number for a wedge. 


It can be shown that Cole’s solution for large negative 
values of goes over exactly into the linearized subsonic 
solution (see appendix A). The reduced-pressure-coefficient 
curve for — 2.02 in figure 7 is so nearly identical for both 
solutions that they cannot be told apart (except that Cole’s 
solution goes to C Vcr at x/c= 1 while the linearized solution 
goes to — oo ) . This is to be expected since the transonic 
perturbation equations are not restricted to transonic 
flow but apply equally well to completely subsonic and 
completely supersonic flow. 3 The transonic equation can 
be written in the form 


Figure 6.— Theoretical reduced local Mach number against reduced free-stream Mach 
number at several chordwise stations on a wedge. Data from references 20 to 22 and 29. 


8 This was pointed out to the author by Dr. Milton Van Dyke of the NACA Ames Aero- 
nautical Laboratory. 
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V 2 “Jdxdx 2 


(19) 


where y is the perturbation potential such that «=f7 + ^> 


v= ~^j' Thus it is clear that for completely subsonic or 

completely supersonic flows the term on the right is negli- 
gibly small but becomes of paramount importance in transonic 
flow. 


CHARACTERISTIC FEATURES OF TRANSONIC FLOW PAST 
WEDGE AND CIRCULAR-ARC SECTIONS 


CHARACTERISTIC FREE-STREAM MACH NUMBERS 

Critical Mach number. — The Mach number at which 
sonic velocity first appears on the wedge is M ga =0 (within 
the inviscid theory) since subsonic flow cannot turn a sharp 
corner. Because of the fact that the boundary layer rounds 
off the corner, and perhaps also because of the spatial 
resolution limitations of the interferometric method, sonic 
velocity was not found there experimentally until 
approximately £<»=-— 0.80 for the wedges. 

The critical Mach number for a half circular-arc airfoil 
followed by a straight section can be obtained approximately 
from linearized subsonic theory. This theory gives the 
surface pressure distribution as 


C v 


— 4 it!c) 
“Wl -MJ 


which vields 


a . 


< 20) 

1.626 (tic) 


VT=m? 


( 21 ) 


at a;/c=0.783 (see appendix B). This equation can also be 
written in transonic similarity form by multiplying both 

( 7 -i-iy/3 

sides by ■ , . . 2 ( 3 — (as shown in the previous section) : 


(t/c ) 2 ' 3 




— 1.626 


( 22 ) 


Now, within the transonic perturbation theory, 


C, = - 2(£-U 

Hence 


(23) 

(24) 


Reynolds numbers used here was fairly thick in comparison 
with the dimensions of the model and thus it “rounded off” 
the shoulder more than would be the case at higher Reynolds 
numbers. (2) The height of the supersonic zone, even for 
an ideal nonviscous flow past thin wedges, appears to be 
quite small until the free-stream Mach number is quite 
close to 1. This is apparent from Cole’s theory and also 
from the argument in reference 25 that the height of shocks 
in the supersonic zone must be of the form 



1/3 A 

c 



(3) The spatial resolution of the interferometric method may 
not have been sufficient to detect very small supersonic 
zones near the shoulder. There is also a large refraction 
error near the shoulder due to the high density gradients 
which tends to obscure details of the flow there. 

Shock-attachment Mach number. — The shock-attachment 
Mach number depends only on the opening angle of the 
profile at the leading edge and can be predicted quite pre- 
cisely by oblique-shock theory. If 6 is the semiopening angle, 
then it can be shown that approximately, for thin profiles, 


M 2 — 1 

00 A 


3 


^ [(t+1)0] 2/3 4 2/3 


(25) 


(see appendix C). If t/c is the thickness ratio of the circular- 
arc section, 0~2(t/c). Hence for the circular-arc profile 



(26) 


Mach number at which sonic velocity appears behind an 
oblique shock. — The Mach number at which sonic velocity 
appears behind an oblique shock M„ s is just slightly higher 
than M mA and again is a function only of the opening angle. 
These values can also be found quite precisely from oblique- 
shock theory and approximately in similarity form can be 
given as 


Ls [(7+ S i)0] 2/3: 
for the wedge (see appendix C) and 

£ =2 


(27) 


(28) 


Equating C Vniu to C Pn one obtains the critical reduced 
Macli number 

{-.,= -0.871 

For the thickness ratio t/c= 0.088 used in these tests, this 
predicts a critical Mach number of 0.834 at z/c=0.783. 
Experimentally, the critical Mach number was found to be 
0.825 and occurred somewhere between z/c= 0.75 and 0.95 
(the pressure distribution was very flat in this range). It is 
interesting to note that the experimental M mcr was higher 
for the wedges than for the circular-arc profile of the same 
thickness ratio. This was probably due to a combination 
of three effects: (1) The boundary layer for the same 


for the circular-arc section. 

CHARACTERISTIC VALUES OF LOCAL MACH NUMBER 

Mach number at leading edge. — The Mach number at the 
leading edge is zero (a stagnation point) for all free-stream 
Mach numbers less than the attachment Mach number. 

Mach number at shoulder of wedge. — The Mach number 
at the shoulder of the wedge just before the turn is always 1. 
This is easily seen in the case of flow with detached shock 
since the only characteristic distance of the finite wedge is 
the distance from the leading edge to the shoulder which 
must determine the shock-detachment distance, and, if the 
sonic point occurred ahead of the shoulder, the shoulder 
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could not influence the shock position. Subsonic flow cannot 
turn a sharp corner so the flow must therefore reach Mach 
number 1 right at the corner. In the case of subsonic free- 
stream flow the argument is not so simple (see reference 22)- 
At the shoulder the flow around the corner is locally a 
centered Prandtl-Meyer fan starting from M = 1 . The Mach 
number just behind the corner is thus determined only by 
the wedge angle and is independent of the free-stream Mach 
number. Behind this point the flow will recompress to the 
free-stream Mach number through a shock or series of shocks, 
for free-stream Mach numbers less than the attachment 
Mach number. The expression for Mach number M P m 
behind an expansion from M= 1 through an angle 6 is 

6 = I tan_1 Vt + T ~ tan "‘ 1 (29) 

Expanding the right-hand side in powers of ^M pm 2 —\ 
(assumed small), the first nonzero term yields 

2(ilW-l) 3/2 

— (30) 

which is in transonic similarity form so 


the Mach number ahead is above 1. This follows from the 
normal-shock relation : 


M x 2 -\ 


i +- 2 J t (m, 2 - 

7+1 


(32) 


1) 


where Mi is the Mach number ahead of the shock and M 2 
is the Mach number behind the shock; so near Mi = l, 

1-M 2 2 «M 1 2 -1 (33) 

or 

1-M 2 «M 1 -J (34) 

Therefore if M O0 = 1 + e, where e is small, the flow past the 
profile is nearly the same as the flow past the profile at 
M oa = l — € since the Mach number behind the central part 
of the detached shock wave is almost exactly 1 — e. It 
follows therefore that the local Mach number distribution 
on the profile surface must have a stationary value at 
M a = l and furthermore vary only slowly in the neighbor- 
hood of M at = l. Mathematically this means 


dM 

dM m 



(35) 


M pm 2 -\ / 3\ 2/3 

[(7+l)0] 2/3 W 


(31) 


PRESSURE DISTRIBUTION ON BODIES MOVING THROUGH 

AN INFINITE FLUID AT SPEEDS NEAR MACH NUMBER 1 

STATIONARY VALUE OF LOCAL MACH NUMBER AT FREE-STREAM MACH 

NUMBER 1 

During the course of these investigations it was found 
that for the wedge and circular-arc sections the local Mach 
number distributions on these sections at very high subsonic 
speeds (above M a , er but below choking Mach number) 
and at very low supersonic speeds (where the detached 
shock wave was a chord length or so ahead of the section) 
were nearly identical. In trying to understand why this 
should be so, the following explanation was derived: (1) At 
low-supersonic speeds the bow shock wave is detached a 
great distance ahead of the profile and a subsonic flow region 
is embedded in the flow field between the shock and the 
sonic line. The part of the shock directly ahead of the 
profile is nearly normal over quite a distance (of course, the 
slope of the shock asymptotically tends to the slope of the 
Mach wave of the free-stream flow at large distances lateral 
to the flow direction). Nagamatsu (reference 30) has pre- 
viously indicated this and points out that the flow past the 
profile should be closely approximated by assuming the 
profile is in a high-speed subsonic flow where the velocity 
distribution at infinity is slightly nonuniform, the minimum 
velocity being directly ahead of the profile and equal to the 
velocity behind the normal shock and then increasing in 
both lateral directions. (2) Now the normal shock near 
Mach number 1 is nearly symmetrical in the sense that the 
Mach number behind the shock is just as much below 1 as 


It should be noticed that this argument is based on two 
assumptions: 4 (1) The detached bow wave moves very far 
ahead of the profile as the flight Mach number decreases 
toward 1. (2) The radius of curvature of the detached bow 

wave at points directly ahead of the profile becomes ex- 
tremely large as the flight Mach number decreases toward 1. 

Examining these assumptions, it would seem that the 
same reasoning should apply to any finite three-dimensional 
body in an infinite fluid traveling at speeds near Mach 
number 1, except that now two radii of curvature at points 
on the detached bow wave ahead of the body must be 
assumed to become large as the flight Mach number de- 
creases toward 1. The detached bow wave is so far away 
from the body at speeds just slightly above Mach number 1 
that the body appears as only a very small object in relation 
to the radii of curvature of the bow wave and, hence, it 
would appear as though the shape and attitude of the body 
could have no appreciable effect in changing the argument 
presented above. 

The reasoning should also apply to an infinite yaw : ed 
cylinder (whose cross section may be finite or, if the angle 
of attack is 0, may extend infinitely far downstream) pro- 
vided that the Mach number considered is the component 
of the Mach number normal to the generators of the cylinder. 

These arguments are for steady-flight speeds. Large 
accelerations through sonic flight speed could conceivably 
modify the phenomenon. Thus it is difficult to judge 
whether or not the available flight-test data confirm the 
concept since nearly all such data come from missile tests 


4 It is believed that these are not actually assumptions but are capable of demonstration 
if one assumes a smooth variation of drag through Mach number 1. 
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that involved large accelerations (or decelerations) through 
sonic flight speeds. The transonic-bump tests of Weaver 
on sweptback wings (reference 1) would seem to support 
these conclusions since they show drag-coefficient maximums 
very near Mach number 1, a necessary consequence of the 
concept for finite three-dimensional bodies and finite, 
unswept, two-dimensional bodies as will now be shown. 

SLOPE OF PRESSURE- AND DRAG-COEFFICIENT CURVES AT Mco = l 

Equation (35) enables one to calculate the slope of the 
pressure- and drag-coefficient curves at Mach number 1 
as follows: 

- p ~ Pa ‘ (36) 



-1 


(37) 


for isentropic flow so 

dC„ 


dM . 


4 

7+1 


2(?p|M a ,=l 

7+1 


(38) 


using 


dM 

dM„ 


= 0 


Now for a two-dimensional body the pressure-drag coefficient 
(based on the chord) is given by the contour integral 


C *=+ <J) C„(i-n)ds 

where 

i unit vector in stream direction 

n unit vector normal to profile pointing outward 

ds element of length along profile contour 

so if the angle of attack is constant and M m is changing 


(39) 


d C d 


dM a 


IX/' 4 2C,\M„-i\ A a 

= 1 7+1 


ds (40) 


But i^)Vnds=0 for a closed contour, so 


dC D 

dM 


2 

7+1 


(41) 


For the front part of a profile (defined as that part ahead 
of the maximum thickness) the usual definition of a drag 
coefficient is 

C Df ,= — - i Cpl-nds (42) 

£ J a 


j: 


where means the counterclockwise line integral from the 


point of maximum thickness on the upper surface to the 
point of maximum thickness on the lower surface; thus 


so 


dCp, 

dM a 

dC Df 


M c 


— A 

= 1 c J. \7+ 1 7+1 / 


dM„ 


4 t 

— COS a- 


7+1 C 


' T + 1 ^ D ^l^oo = l 


( 43 ) 


where 

t maximum thickness of profile 
a angle of attack of profile 

Similarly the drag coefficient for the rear part is 


so 


dC Dt 


dM a 


Cd= — ~ f Cpi-nds 
c Jb 


-4 t 2 n 


(44) 


(45) 


For the tests on wedge and circular-arc sections followed 
by straight sections the concept of drag coefficient of the 
front part of the section will often be used. 

For bodies of revolution (which include spheres, cone- 
cylinders, etc.) the pressure-drag coefficient (based on maxi- 
mum cross-sectional area) at zero angle of attack is 


Cr 


= L+ P< *(fl) 


(46) 


where 

R maximum radius of body 
l length of body 
x distance from nose along axis 
Therefore 


dC D 


dM a 


— 2 C d \. 


7+ 1 


as before in the two-dimensional case. However, for front 
and back drag coefficients 




SO 


dC d. 


dM 

and similarly 

dC D 


dM n 


4 2 

± r n 

"7+1 7+1 D Ha/oo = i 


4 2 

- — f! I 

7+1 7+I 


(4 7) 


(48a) 


(48b) 


and these differ from two-dimensional values obtained above 
in equations (43) and (45) by not involving the fineness 
ratio of the body (this is of course due to the different 
reference areas for drag coefficients). 

For the general finite three-dimensional body the pressure- 
drag coefficient is given by 
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Gd=~\ ^C v {%-rt)da (49) 

where A is some reference area of the body and S is the 
surface of the body. It follows as it did previously that 


Mach number distribution near M a3 = 1 for airfoils at an 
angle of attack is also shown clearly in figures 8, 9, and 10 
of reference 32. 

ON COMPARING THEORY AND EXPERIMENT 


dC D 
dM„ 

SLOPE OF DRAG-COEFFICIENT CURVE AT Moo = 1 IN TRANSONIC SIMILARITY 
PARAMETERS FOR TWO-DIMENSIONAL FLOWS 


7+1 


C D I 


Within the transonic approximation 

C v =- 2(£-£j 

so 



(50) 

( 51 ) 


xt dM 
Now -j-ttt- 
dM m 


Now 


so 




= 0 implies that 4^-1 =0; hence 

= 1 «£j*oo=0 

dC v \ 


d £ 00 |foo=0 


= 2 


C p i'hd 


(9 


dC D 


= 0 


£=o=0 


Similarly it is easy to show that 

dC Dp 


and 


dL 

d_C DR 

dL . 


f»=o 


= 2 


£»=0 


= — 2 


(52) 

(53) 

(54) 

(55) 

(56) 


OTHER DATA SHOWING SLOW VARIATION OF LOCAL MACH NUMBER NEAR 

Moo = l 


As mentioned previously, Maccoll in 1946 had already 
proposed the slow variation of local Mach number near 
M m = 1 on “bodies having distinct corners . ” It appears 
that this latter restriction is not necessary. MaccolPs pro- 
posal was based on rather slim evidence and it is believed 
that here, on the basis of the argument presented concerning 
the normal shock, the principle is explained more convinc- 
ingly. Also the experimental evidence given here and by 
Drougge (reference 18), Bleakney and Griffith (personal 
communication), Weaver (reference 1), and by some NACA 
reports tends to bear out the conclusions of slow variation of 
local Mach number on bodies near M m = 1. 

This fact is sometimes slightly obscured in the NACA 
reports because pressure coefficient was plotted instead of 
plpo or local Mach number. However, constant Mach 
number lines were sometimes drawn in these plots and there 
the evidence shows up strongly (see, e. g., reference 31, figs. 
7 to 11, pp. 36 and 37). The relative constancy of local 


In references 21 and 25 discussions were presented on the 
philosophy of comparing experiments with approximate 
theories, and these discussions will not be repeated here > 
except to mention that in some of the theoretical curves 
presented here the values have been shown with a certain 
spread which results from using a pressure coefficient equal 

to —2 s- 2 or —2 — ~ (the former value is the one that 

a U 

fits into transonic similarity theory; the latter value is the 
one more commonly used in perturbation analysis). 

In connection with the idea presented in reference 25 of 
extrapolating experimental data to zero thickness in order 
to compare with results from transonic perturbation analyses, 
it is interesting to note that the characteristic Mach numbers 
mentioned in the section “Characteristic Features of Tran- 
sonic Flow past Wedge and Circular-Arc Sections” can be 
presented in powers of the thickness of the wedge (or equiv- 
alently in powers of the wedge angle), the first term of 
which gives the transonic similarity expression; two of these 
values are 

M 2 — 1 
*” s= [(y+1)0] 2/8 

= 2 1/3 [l + 19 ^~ 2 (^i)‘ /3 fl^+^fl 4 ' 3 )] (57) 

y M pM 2 1 

*™ - [(7 + 1)0] 2/s 

=(lH'+^[^T # “’ +0< ' n i (58) 

(See appendixes D and E.) In transonic perturbation 
theory the terms in 0 on the right-hand side are neglected. 
This can lead to fairly large errors for even moderately 
large values of 0 since the approach to 0=0 is nonlinear and 

< 59 > 

Judging from this one might expect that quantitative 
agreement of transonic perturbation analyses with experiment 
would not be so good. However, in comparing two similar 
shapes with only slightly different thickness ratios by 
transonic similarity considerations one would expect fairly 
good agreement. 


EXPERIMENTAL RESULTS 
FLOW FIELD NEAR 10° WEDGE 

Figure 9 shows interferograms of the flow past the 10° 
semiangle wedge for 14 Mach numbers from 0.700 to 0.892 
and 1.207 to 1.465 (the interferograms for the 4j4° and 
7 %° wedges were very similar and hence are not shown here). 
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Figure 9.— Interferograms of flow past 10° semiangle wedge for various Mach numbers. 



(c) Aia =0.852; p m /p 0 = 0.712. 


Figure 9.— Continued. 


(d) Af a, =0.892; Pco/po=0.691, 
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(f) 


(e) Moo = 1-207; Pco / Po =0.528. 


(f) Moo =1-240; p oo /p o =0.511. 


Figure 9.— Continued. 




(h) 


tg) Me, =1.278; Poo / Po =0.493. 


(h) 3 /oq =1.315; p co /po=0.476. 


Figure 9.- Continued. 


243695-53 3 


16 


REPORT 1094 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



(i) =1.356; p cx> /po=0.457. 


Figure 9.— Continued. 


0) A/co =1.375; p oo /p o =0.449. 



(k) 





(k) =1.391; p„/p.=0.441. 


Figure 9.— Continued. 


(/) Mo, = 1.411; Poo/p.=0.432. 
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Notice that the lines of constant density in t lie subsonic-flow 
interferograms are roughly elliptical in shape as predicted 
by the theory (see appendixes A and B). A supersonic 
flow region was first detected between = 0.700 and 0.794 
(the sonic line is shown as a dashed line in the figures) and 
a shock emanating from the corner appears in the super- 
sonic zone at = 0.794. As the Mach number was 
increased, this zone grew larger and a shock appeared at 
the rear of it, while the shock emanating from the corner 
weakened and disappeared. This rearward shock was of 
the typical X type associated with a laminar boundary 
layer, and the interferograms clearly indicated the separation 
of the boundary layer ahead of this shock. The similarity 
between the flow field at M m = 0.892 and at M^l.207 
(figs. 9 (d) and 9 (e)) is striking; the base of the rearward 
shock has moved quite far back on the wake of the blunt 
trailing edge at M m = 1.207 but in the vicinity of and ahead 
of the sonic line the two fields are nearly identical except 
for the detached shock wave which appears about lK chord 
lengths ahead of the wedge at M m = 1.207. As the Mach 
number was increased above 1.207, the detached shock 
moved in closer to the leading edge and finally “attached’ 7 ' 
at a Mach number quite close to the theoretical attachment 
Mach number of M m = lA\8. Notice that the process of 
attachment is very continuous. The effect of the boundary 
layer is quite noticeable in the last few interferograms: 
This can be roughly accounted for by considering the 


boundary layer to clitage the shape of the body by its 
displacement thickness and then considering a nonviscous 
flow past this revised shape. On the wedge the boundary 
layer will not grow so rapidly as on a flat plate because of the 
favorable pressure gradient and, in fact, the effect of the 
strong expansion around the corner is known to cause an 
almost complete collapse of the boundary layer there. As 
the bow shock wave gets close to attachment, the velocities 
in the subsonic region behind it are getting very close to 
sonic velocity and hence the flow in this region is very 
sensitive to any slight curvature of the “revised shape” 
of the wedge. This accounts for the shift of the base of the 
sonic line forward to the leading edge as the shock ap- 
proaches attachment. The nonviscous theory would in- 
dicate that the sonic line would always begin at the corner 
and, at a Mach number just slightly above the shock- 
attachment Mach number, the whole subsonic region would 
become sonic; then, with increasing Mach number, the flow 
behind the shock would be completely supersonic. As 
observed, the boundary-layer effect is to make the wedge 
have a curved surface and the sonic line actually moves 
slowly from the corner to the nose. Even with an attached 
shock wave at M m = lA6o the flow behind the shock is not 
quite uniform (as nonviscous theory would indicate it 
should be) because of the effective curved surface caused 
by the boundary layer. 
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LOCAL MACH NUMBER DISTRIBUTIONS ON THREE THIN WEDGES 

Figure 10 shows the variation of local Mach number 
distribution on the surfaces of the 4%° , 7%°, and 10° semiangle 
wedges with free-stream Mach number. This should be 
compared with figure 5 which shows the corresponding 
theoretical curves in terms of the transonic similarity param- 
eters. The general behavior of the theoretical and ex- 
perimental curves is quite definitely in good agreement. 
Particularly noteworthy is the slow variation of the local 
Mach number distribution near free-stream Mach number 1. 

PRESSURE-COEFFICIENT DISTRIBUTIONS ON THREE THIN WEDGES 

The slow variation of the Mach number distribution in 
the range near M m = 1 is obscured when the results are 
plotted in terms of pressure coefficient, since the pressure 
coefficient changes a great deal if local Mach number is con- 
stant while the free-stream Mach number changes. A better 
parameter for presenting transonic pressure distributions 
would be plpo (plpo in case of a detached shock). Typical 
C p distributions are shown in figure 11 for the 7%° wedge 
(the results for the 10° and 4 %° wedges were very similar 
and hence they are not presented). The points shown were 
where the fringes intersected the body in the interferograms. 



0 .2 .4 .6 .8 1.0 1.2 

x/c 

(a) 10° semiangle wedge. 

Figure 10.— Local Mach number against x/c for increasing free-stream Mach number. 



(b) 7h° semiangle wedge. 
Figure 10— Continued. 



(c) 4 Vf semiangle wedge. 
Figure 10.— Concluded. 
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0 .2 .4 .6 .8 1.0 


x/c 

(a) Subsonic. 

Figure 11.— Pressure distributions for 7^° semiangle wedge. 


Since for a wedge the drag coefficient is proportional to the 
average C P) the drag rise is evident in the subsonic distribu- 
tions as the point where C p =0 moves rearward with increasing 
free-stream Mach number. Linearized subsonic theory 
(which predicts C D = 0) locates the C p = 0 point at x/c=5 0 
percent. Figure 7 shows theoretical reduced C p distributions 
at various reduced free-stream Mach numbers. Again the 
qualitative agreement of these curves with experiment is 
evident. 


SHOCK-DETACHMENT DISTANCE FOR THREE THIN WEDGES 


Figure 12 shows the shock-detachment distance against 
reduced free-stream Mach number for the three thin wedges 
and includes the theoretical values from reference 21. Here 
Vincenti and Wagoner’s values for have been multiplied 


by 


^A\d=7H° 

^ °°i4 1 0 = 0° 


in order to make the transonic perturbation 


value of detachment reduced Mach number agree with the 
value from oblique-shock theory for the wedge: 5 The 
reason for this was discussed in the section “On Comparing 
Theory and Experiment,” namely, the difficulty of comparing 




(b) Supersonic. 
Figure 11. — Concluded. 


transonic perturbation theory quantitatively with experi- 
ment. Notice how rapidly the shock wave moves away 
from the wedge as the Mach number is decreased toward 1 . 

DRAG-COEFFICIENT VARIATION WITH MACH NUMBER FOR THREE THIN 

WEDGES 

It was shown in reference 25 that the viscous effects on 
the wedge tend to compensate each other at the leading 
edge and the shoulder so that the over-all pressure drag is 
nearly the same as if the flow were inviscid. Thus it would 
be expected that the pressure-drag coefficients obtained by 
integrating the experimental pressure distributions would 
check the inviscid transonic perturbation theory. The 
reduced drag coefficient used here was 

^-JX!) (6o > 

which is, in essence, the reduced drag coefficient of the 
upper (or lower) half wedge. This was done since the wedge 
model was regarded as the front half of a double-w^edge 
profile and lienee the value given here is the part of the re- 


5 In terms of Mach number, for the 7\h° wedge the shock theory predicts attachment at 
Mm ,4 = 1-33 (Z<x> A = 1.68), while the transonic perturbation theory predicts Moo x = 1.25 (See* 
= 1.19). 
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Figure 12.— Shock-detachment distance against reduced Mach number for a wedge. 


duced drag coefficient contributed by the front half of such 
an airfoil ( C Dp as in equation (42)), based on the chord of the 
double-wedge profile, which would be twice the chord of the 


model used here. Of course, this viewpoint is valid only for 
supersonic free-stream Mach numbers. 

Figure 13 shows the reduced drag coefficients for the three 
thin wedges plotted against reduced Mach number. It is 



semiangle wedge] - 
7^' semiongle wedge Experim " ent 
4 semiangle wedge J i . 

10° semiangle wedge) ^hock 
7 ¥ semiangle "edge l extension 
4-^ semiangle wedgej theory 


Theoretical slope at £ =0 

j i > 


3.0 


Figure 13.— Experimental reduced drag coefficient against reduced Mach number on a wedge. 


seen that the results give nearly a universal curve, which 
they should if the transonic similarity law is true, but that 
there are systematic variations with wedge angle. This is 
to be expected based on the discussion of the section “On 
Comparing Theory and Experiment.” The vertical lines 
through the experimental points indicate estimated accuracy 
of the data. This figure should be compared with figure 8, 
the theoretical reduced-drag-coefficient variation with re- 
duced Mach number. It is obvious that the qualitative 



Figure 14.— Drag coefficient against Mach number for 4^°, 7 \^°, and 10° semiangle wedges. Comparison of theory with experiment. 
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agreement of theory and experiment is good. In figure 14 
the theory and experiment are compared directly for the 
three thin wedges. Here the theoretical drag coefficients 
are shown with a vertical spread, the upper values for 
1 corresponding to the use of the pressure coefficient 

C n = — x — — and the lower values, to the use of the pres- 

p a * 

sure coefficient C p = — The situation is vice 

versa for From this figure it is evident that the 

transonic perturbation theory gives a good approximation 
to experiment. 


FLOW FIELD AT =1.44 FOR THE 26.6° WEDGE 

Figure 15 shows the experimental and theoretical constant- 
velocity lines in the subsonic region behind the detached 
shock wave for a 26.57° semiangle wedge at M 00 = 1.44. 



The theoretical analysis was made from relaxation calcu- 
lations by Drebinger (reference 19) who used the flow equa- 
tions with entropy variation behind the shock taken into 
account. The experimental constant-velocity lines were 
determined from the isopycnic lines of the interferogram by 
taking into account the lateral stagnation-pressure gradient 
behind the curved shock. The isopycnic lines near this 


strong shock wave were probably slightly in error because of 
the “smearing out” of the pressure discontinuity across the 
shock in the side-wall boundary layers. It is seen that the 
agreement between theory and experiment on detachment 
distance and constant-velocity contours near the wedge is 
good . 

Figure 16 shows the surface pressure distribution from 
reference 17 and the present experiments. Again it is seen 
that the agreement is good. 



x/c 

Figure 16.— Pressure distribution on a 26.6° semiangle wedge at = 1.440. 

FLOW FIELD NEAR THE 8.8-PERCENT CIRCULAR-ARC SECTION 

Figure 17 shows interferograms of the flow past the 
8. 8-percent circular-arc section for 14 Mach numbers from 
0.718 to 0.936 and 1.11 to 1.500. 

Supersonic velocity first occurred at A/*, = 0.825 (see the 
section “Critical Mach number”) and in figures 17 (c) and 
17 (d) a nearly symmetric supersonic zone is shown at 
M m = 0.848. No shock waves were apparent in this zone, 
although a sensitive schlieren apparatus might have shown 
some weak shocks there. At M m = 0.890 the supersonic 
zone has grown rapidly and now terminates in the X-shock 
configuration. Further increase of the Mach number to 
Mo^O.935 (figs. 17 (e) and 17 (f)) shows the supersonic 
zone increasing laterally and the terminating shock moving 
rearward into the wake of the body. Figures 17 (e) and 
17 (f) also show the density distribution at (the 

detached shock wave was just out of the field of view of the 
interferometer) and it is interesting to note the similarity 
between the flow field at M m = 0.935 and 1.11. It 

would appear as though the shock terminating the super- 
sonic zone at M a) = 0.935 had moved rearward to form the 
trailing-edge shock (which is actually in the wake here 
because of the blunt trailing edge) and the supersonic zone 
had grown laterally until the sonic line joined with the 
detached shock far away from the body at M m = 1, thus 
causing an embedded subsonic zone in the supersonic flow 
with further increase in Mach number. 
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.795-. 



(o) 



(b) 


(a) =0.718; Poo / Po =0.783. 


(b) JV/qj =0.819; Poo / Po =0.730. 


Figure 17.— Interferograms of flow past 8.8-percent circular-arc section for various Mach numbers. 



(c) A/oo =0.848; Poo / Po =0.715. 


Figure 17.— Continued. 


(d) Moo =0.890; Poo / Po =0.696. 
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(e) M oo =0.935; Poo / P ,=0.668. 


Figure 17.— Continued. 


(f) Moo = 1.11; Pm/p o=0.577. 




(g) =1.160; Poo / P o=0.551. 


Figure 17.— Continued. 


(h) M m =1.200; Poo / Po =0.531. 
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Figure 17.— Continued. 



Local Mach number, 
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(m) M m =1.450; p oo /p ( ,=0.416. (n) =1.500; p «//>«= 0.395. 

Figure 17.— Concluded. 



Figure 18.— Local Mach number against x/c for increasing free-stream Mach number. 8.8- 
percent circular-arc section. 


With further increase of Mach number above M m = 1.11, 
figures 17 (g) to 17 (n) show that the detached shock again 
approached the leading edge and the embedded subsonic 
zone decreased in size until finally the shock “attached” 
somewhere between M 00 = 1.400 and M oo = lA50 (the theo- 
retical value being M mA = 1.423). 

LOCAL MACH NUMBER DISTRIBUTIONS ON 8.8-PERCENT CIRCULAR-ARC 

SECTION 

Figure 18 shows the local Mach number distributions for 
the 8. 8-percent circular-arc section as obtained from the 
experiments at various free-stream Mach numbers. Again 
it is apparent that the variation of local Mach number 
distribution near AI 00 = 1 is very slow and, indeed, the dis- 
tribution for M m = 1 could be interpolated from this figure 
with good accuracy. 

Figure 19 is a cross plot of the data of figure 18 except 



Figure 19. — Variation of reduced local Mach number distributions with reduced free-stream 
Mach number. 8.8-percent circular-arc section. 
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that here the data are given in transonic similarity param- 
eters. This figure shows contours of constant reduced 
local Mach number on a plot of reduced free-stream Mach 
number against chordwise position. The dashed lines repre- 
sent subsonic local Mach numbers; the solid lines, supersonic 
local Mach numbers. Note again the slow variation of 
local Mach number distribution with free-stream Mach 
number near sonic velocity. 

PRESSURE-COEFFICIENT DISTRIBUTIONS ON 8.8-PERCENT CIRCULAR-ARC 

SECTION 

Figure 20 shows the pressure-coefficient distributions on 
the 8.8-percent circular-arc section for various free-stream 
Mach numbers. The points shown are where the fringes 
intersected the body in the interferogram. Again the 
presentation in this manner obscures the interesting fact 
observed in figure 18. 

DRAG-COEFFICIENT VARIATION WITH MACH NUMBER FOR 8.8-PERCENT 
CIRCULAR-ARC SECTION 

Figure 21 shows the experimental determination of the 
drag coefficient of the front part of the 8. 8-percent circular- 
arc section. This again is of the nature of a fore drag coef- 
ficient and, as shown in equation (43), it should have a 



(a) Subsonic. 

Figure 20.— Pressure distributions on an 8.8-percent circular-arc section. 



x/C 

(b) Supersonic. 

Figure 20.— Concluded. 

4 t 2 

positive slope equal to — ^ C D | A / 00 =i at M 00 = 1 ; this 

y \ i c 7t i 

is how the subsonic data have been joined with the super- 
sonic data. The vertical lines through the experimental 
points again indicate estimated accuracy of the data. For 
the case of an attached shock the pressure distribution can 
be calculated using characteristics theory and the shock 
polar; however, a close approximation is obtained by con- 
sidering the flow behind the shock wave to be Prandtl- 
Meyer flow. (This yields, approximately, parabolic-shaped 
bow and trailing-edge shock waves; see reference 33.) From 
this pressure distribution the drag was calculated and is 
shown in figure 21. Taking into account the “reflected” 
characteristics from the shock wave would give more com- 
pression and increase the drag coefficient so that it would 
agree better with the experimental values at M CO = 1A50 and 
1.500 shown in figure 21. 

Note that the tests were made at low enough supersonic 
speeds to get definitely below the drag-coefficient maximum 
at M m ~ 1.20. 
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LOCAL MACH NUMBER DISTRIBUTIONS ON A 12-PERCENT BICONVEX 
CIRCULAR-ARC AIRFOIL 

Figure 22 shows local Mach number distributions from 
reference 8 for high -subsonic-speed flow over a 12-percent 
biconvex circular-arc airfoil (with turbulent boundary 
layer). The data for the 8.8-percent circular-arc section at 
two supersonic speeds have been scaled according to the 
transonic similarity laws to the 12-percent case and are 
shown for the front half of the 12-percent airfoil in figure 22. 
The back half for these two cases has been faired in using a 
Prandtl-Meyer expansion which should be approximately 
correct (a more accurate determination could have been 
made using characteristics theory and the shock polar). 
At M 00 = 1.58, the theory indicates that the shock is attached 
with sonic speed just behind the shock on the leading edge, 
so that the distribution can be obtained by standard methods 
mentioned above; again the Prandtl-Meyer expansion ap- 
proximation was used for the distribution at M 00 = 1.58 in 
figure 22. 

The behavior of the Mach number distributions is similar 
to that of the distributions shown previously, except in this 
case the movement of the shock terminating the local super- 
sonic zone is shown. Apparently little change in local Mach 
number distribution occurs between M m = 0.936 and 
M m = 1.29. 

DRAG-COEFFICIENT VARIATION WITH MACH NUMBER FOR A 12-PERCENT 
BICONVEX CIRCULAR-ARC AIRFOIL 

The data of figure 23 were converted to pressures which 
were integrated to give the pressure-drag coefficient for the 
various free-stream Mach numbers. In addition the drags 
of the front and back halves are shown separately. The 
drag-coefficient variation between M oo = 0.96 and M O0 = 1.20 
was based on constant local Mach number distribution at 
values interpolated between the curves for = 0.936 and 



Figure 22.— Local Mach number against x /c for increasing free-stream Mach number. 12* 
percent biconvex circular-arc airfoil. 
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M m = 1.29. The data were faired into the curves for attached 
shock wave calculated on the Prandtl-Meyer expansion 
basis. It is seen that the fore drag coefficient has a maximum 
after = \ while the drag coefficient of the rear part has a 
maximum before M m — 1. The over-all airfoil has a maxi- 
mum drag coefficient just before M O0 = l in order for the 
curve to have the slightly negative slope at M m = 1 given 
by equation (41). 

CONCLUSIONS 

An experimental investigation of transonic flow past two- 
dimensional wedge and circular-arc sections was made using 
a Mach-Zehnder interferometer. The conclusions may be 
stated as follows: 

1. The transonic similarity theory of Von Karman and 
Guderley was checked and found to be in good agreement 
with experiment for thin wedge profiles near a free-stream 
Mach number of 1. 

2. The results of theoretical calculations, using transonic 
perturbation theory, made by Guderley and Yoshihara, 
Vincent i and Wagoner, and Cole for a wedge in transonic 
flow were checked experimentally at high -subsonic and low- 
supersonic speeds for three wedges of different angles and 
were found to be in good agreement with experiment. 

3. The flow field and the surface pressure distribution for 
a 26.6° semiangle wedge at a free-stream Mach number of 
1.44 were obtained experimentally and were found to be in 
excellent agreement with the theoretical calculations of this 
flow made by Drebinger. 


4. The pressure distributions and drag coefficients for an 
8.8-percent circular-arc section followed by a straight section 
and for a 12-percent biconvex circular-arc airfoil were pre- 
sented completely through the transonic range. It was 
shown that some difficulty arises in comparing two-dimen- 
sional transonic perturbation theory with experiment, since 
this theory neglects thickness-ratio terms of order ( t/c ) 2/ 3 
and higher; for even moderate thickness ratios this will 
cause noticeable deviations from more exact theory. 

5. It was shown from some physical arguments that the 
local Mach number distribution on bodies traveling through 
an infinite fluid has a stationary value at free-stream Mach 
number 1. This was verified experimentally for the case of 
two-dimensional flow. It was shown that this concept implies 
a drag-coefficient maximum just below free-stream Mach 
number 1 for all bodies in steady flight. This fact can be 
used to obtain the variation of local Mach number distribu- 
tion on bodies completely through the transonic range of 
velocities from wind-tunnel tests, provided small models 
are used so that tests can be carried well above critical Mach 
number and to low enough supersonic Mach numbers so 
that the bow shock wave is detached a chord length or so. 


California Institute of Technology, 
Pasadena, Calif., June 1 , 1951. 
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ASYMPTOTIC REPRESENTATION OF COLE S SOLUTION FOR LARGE NEGATIVE VALUES OF REDUCED FREE-STREAM 

MACH NUMBER 


Cole’s solution for the high-subsonic-velocity flow past a 
thin wedge (reference 22) is given as follows (in Cole’s 
notation) : 




sinh \(v 0 —v) 
sinh \ v 0 


J - i/3(Xz\)\d\ 


(61) 


x(Z } V]Zi) 


= l- Zl m 2 m Vo j“ o) J 2/3 (\z)J_ m (\ Zl )\d\ 

(62( 

where the center line of the wedge is at y= 0; the leading 
edge, at x=0; the shoulder, at x—1) and 

13/2 


Z\ 




tf o =(7+l)0 


2 \ 

I 


(63) 


and the other notation is the same as that in the present 
paper. 


Using the standard methods of partial-fraction expansion, 
one may write 

sinh \{v 0 — v) 0 ” . / v\ \ 2 v 0 2 , aA , 

7. - 1 ~ S 8111 (r V.) <64) 


sinh X 
cosh \(v 0 — v) 1 


sinh \ v 0 


\v 


- 2 X) cos 


O 71 = 1 


(ut— ^r-2 2 2 (66) 

V vJ\ 2 v 0 2 +nV 


( 66 ) 


Substituting these into the integrals above and making use 
of the integrals 

f \ 2 JT ~2 J-i/s(Xp)J-i /z (\y)d\ = 

(-ot 2 I_ l/3 (ay)K- m (ap)] 0 > 7 >O 
\ — oc 2 Ki, z (ay)I- m (aP)] 7>/3>0 

and 

F \2 I 2^~ 1 /^P)^2/z(Xy )d\ = 

Jo A -f-a 

( — aI 2 ,%{ay)K-iiz{afi)) /3>7>0 
\otK 2 iz(ay)I-i,z(aP); 7 >/ 3>0 
equations (61) and (62) can be written as 


(67) ( 


y=v 0 m y= 
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Making use of the asymptotic formulas 
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I’'( z )~~fiT = Z e *+ ■ ■ • as 
y Ztz 


K,(z)~Jf z e-’+ . . . 

and the simple summations 

■ sj — \ —nn • oil! /I U 

Zj « sin mrx= 7i r » 

,7Ti 2 cosh a — 2 


as z —- > <» 


sin it a; 


COS ITJ' 


(70) 

(71) 

(72) 


cos wx—e 


Ec an cos nirx^-r. , 0 

Jfrt 2 cosh a — 2 cos irx 


(69) 7 


(73) 


one can write equations (68) and (69) for large values of z 
and Zi as 


-Kintio' 


1/6 


V 

sin 7 r — 

Vo 


, 7r(2 — 3i) 2/ 

cosh — cos 7r 


; 2 , Zi->oo 


( 74 ) 


9 Equation (67) is from reference 34; equation (66) is obtained by differentiation of equation (67) with respect to 7. 

7 Figures 5, 6, and 7 were calculated from this equation for p=0, for the cases £„ ^0. 
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X ~ < 


1 + 


, , N 1/6 COS 7T £*<*-», >/«» 

1 /2\ 1/f v 0 


m 


x 7 r(z — Zi) V 

cosh — - — COS 7 r 


o) > 00 


Vo 


Vo 


(75) 


i/ \i/6 costt 

_ 1 /_ 2 _\ 1/6 v 

2 \zj 


, 7r(2 — 2i) V 

cosh — - — cos 7 r 


0 ; z>2i 


V 0 V 0 

Eliminate v between equations (74) and ( 75 ) for Zi>z— >o° 

(z/zO^V 


( 


y 


v ( , y 

(_i_) "^2 sinh [ir (3 — Zi)/v 0 ] ( , 

/ I \ / / \1 /ft / i 


) (z/z i) 1/6 r) (2/2i) 1/6 f 

(2(3 Zi/2) 1/3 sinh [tt(z — z i)/ 0 «j) ( 2 sinh [tt(z — Zi)/»J ) 

(76) 

and for z>Zx— >oo simply replace z by z — 1 and z— Zi by 
Zi— z in equation (76). Thus the lines of constant Mach 
number are ellipses with centers on y=0, with ratio of semi- 
axes equal to 

(77) 

which is precisely the solution given by the linearized 
subsonic theory (see appendix B). 

Now, in the notation of the present paper, 

T(z ~ Zl) =~ [(- *) 3/2 -(- U 3 ' 2 ] 

Vo o 

and since £— is small on the wedge and since 

c p =- 2 (a-u 


one can write 




Hence, 




so, approximately, 


ir( z — 2,)_ IT, y MltTl _ TT Cp^l—MJ , s 

c > — 2 e 2 ’ 2l ~*°° (78) 

Similarly, for large values of 2 and Zi it follows that 




l;z,Zj— >00 for z — Zi small 


(79) 


Substituting equations (78) and (79) into equation (16), one 
gets the exact linearized subsonic solution for constant- 
velocity lines (see appendix B). Therefore on the wedge 
(y= 0), from equations (76), (78), and (79), one has approxi- 
mately for large values of z and 




,V^c p 

e — e 


Solving this for C v , 


c , V -*.= -§ lo & V 5 


c v =- 


-20 
TT \1 — 


l0£ 


x-l 


(80) 


which is precisely the linearized subsonic solution for flow 
past a wedge (see appendix B). Thus Cole’s solution far 
away from M ao = l tends exactly to the linearized subsonic 
solution. 


APPENDIX B 

LINEARIZED SUBSONIC AND SUPERSONIC FLOW PAST WEDGE AND CIRCULAR-ARC SECTIONS 


LINEARIZED SUBSONIC FLOW PAST A WEDGE 

Let the wedge center line be on y= 0 , with the leading 
edge at £=0 and the shoulder at x=c. Then the incom- 
pressible-flow problem is to find an analytic function u—iv 
such that s ;=0 on y= 0 except for 0 <z<Cc where v=Ud 
and u—iv= 0 at infinitv. Such a function is 


UO , z/c 
u — iv = — log* t — 7 
7 T & (z/c )— 1 


(81) 


where z=x-\-iy and U is the free-stream velocity. Thus on 

y= 0, 

-2 u —26, x/c /oox 

TT = -ir log ‘T=w5) (82) 


Cp 


Using the Prandtl-Glauert transformation, for linearized 
subsonic flow 


Op 


Cp 


VI -MJ ir\ 1 — A/ 
or in transonic similarity notation 


20 x/c 

log 


1 — (x/c) 


(83) 


C r = j= loge - X/C - 
£» 1 — (x/c) 


(84) 


For the incompressible case the lines of constant pressure 
in the fluid will be where 


z/c 


:= Constant 


(z/c)- 1 

but these are circles v T ith centers at 


and radii 


-rC, 

1 — e 6 


y = o 


(85) 


— * Cr 

o 2 0 


-*C P 


1—e 


In the Prandtl-Glauert transformation the y distance is 
stretched by the factor Vl — as is the pressure coeffl- 


[ 


y 


TRANSONIC FLOW PAST TWO-DIMENSIONAL WEDGE AND CIRCULAR-ARC SECTIONS 

On the wedge (y = 0, then 

so the linearized subsonic solution is 


cient so the lines of constant pressure (and hence density) 
are ellipses with ratio of axes equal to -Jl—MJ 2 given by 
the equation 


7 


+ 


l/2Vl ~MJ sinh 


2 sinh (irC^y 1 — M+/2 O) 
l/sinh (irC^l-Mj/26) 


= 1 ( 86 ) 


LINEARIZED SUPERSONIC FLOW PAST A WEDGE 

From the Ackeret theory the pressure coefficient in super- 
sonic flow is proportional to the slope and for the wedge 
yields simply 

Cp = -^== (87) 


or 


p 4L 


1 


( 88 ) 


LINEARIZED SUBSONIC FLOW PAST A CIRCULAR-ARC SECTION 

For the circular-arc section, the slope of the surface 
varies almost linearly with distance from the zero-slope 
point along the axis of the profile. For the section shown 
in figure 1 then, with the center line on y= 0, the leading 
edge at a;=0, and the zero-slope point at x=c, the incom- 
pressible-flow problem is again to find an analytic function 
u—iv such that on y=0, v=0 except for 0<£<c where 

r — ( l — — j where t is the half thickness at x=c and 
c\ cj 


--2U 


u —{ v — o at infinity. Such a function is 

v![(f— O' 08 - (7/oRi- 1 ] < 89) 
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(90) 


c ’=-vr 


4 t_ 

7T C 


or 




(91) 


(92) 


The minimum C v is obtained by differentiation, and one 
finds that 

dC v 


at the point where log* 


d(x/c) 

x/c 


-0 


1 


1— (x/c) x/c 
solution of this transcendental equation is 


and numerically the 


which gives 


x /c= 0.783 
— 1.626 (t/c) 


Cpmtn yl-MJ 


(93) 

(94) 


LINEARIZED SUPERSONIC FLOW PAST A CIRCULAR-ARC SECTION 

The result here is again simple from the Ackeret theory: 

4 m 


Cp 


or 


C„ 


= Vt 


h('~) 

fcH) 


(95) 

(96) 


APPENDIX C 

TRANSONIC SHOCK POLAR 


The equation of the shock polar in the hodograph plane is 


u 

T*“ 


V 2 = (U — U) 2 —jr 


U 


u 


u 

1 a*' 


u 

T* 


(97) 


wffiere U is the velocity ahead of the shock and u and v are 
velocity components behind the shock parallel and perpen- 
dicular to the direction of U, respectively. Making the 
transonic approximation in this equation, let 


(98) 


Substituting into equation (97), neglecting higher powers of 
the perturbation velocities, one obtains 





(u„-«) 2 (J+^) 


(99) 


Letting 


one then has 


u' = ( 7 + 1 ) 
t»'=(r+l) 
2{vJ={uJ-uJ{u'+uJ) 


n 

pf 


( 100 ) 


( 101 ) 


The wedge angle for detachment of the shock will now 
be given by the maximum value of v r . This is easily seen 


to occur at u'= —~r u a 

o 

within the transonic approximation 

v'=(y+l)6 

1 


giving iw'=gV3 uj n . Since 


on the wedge and 
this implies that 

(7+1)0* 


• 3 3/ 2 


(A/„ 2 -l) 


, 3/2 


( 102 ) 

(103) 


or, viewed in another light, this implies that the reduced 
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attachment Mach number is 

l 




A [(t+1)0] 2/3 4 2/3 


1.19 


(104) 


Similarly, the wedge angle for obtaining exactly sonic ve- 
locity behind the shock is given by the value of v' where 

u'= 0 . This is v'=—j= (uj) 312 . Again using equation ( 102 ), 

y2 


this can be written 


1 


( 7 + 1 ) 0 . 1) 3 ' 2 

V2 


(105) 


and, viewed in another light, this implies that the reduced 
Mach number for which sonic velocity is obtained behind 
the shock on a wedge is 


L 


M„ s -l 


[(7+1)0] 


2/3 


: 2 1/3 = 1.26 


(106) 


APPENDIX D 


VARIATION OF REDUCED MACH NUMBER AT WHICH SONIC VELOCITY OCCURS BEHIND AN OBLIQUE SHOCK WITH 

FLOW-DEFLECTION ANGLE 


The oblique-shock relations can be written 
7—1 


1- 


Mi 2 


M 2 cos 2 0 


7 M 2 sin 2 0 — * 1 + 7 — 1 M/ sin 2 0 


tan 0=2 cot 0 


2 1 2 
M\ 2 sin 2 0—1 
Afi 2 (7+cos 20) + 2 


(107) 

(108) 


where 


Mi Mach number ahead of shock 
M 2 Mach number behind shock 
0 shock-wave angle 
0 flow deflection angle 

For M 2 = 1, eliminating 0 between these two equations yields 


tan 0 = 


/(MO 


-yMi 2 -f(Mi)~ 


M, 2 -l- 7 9 + 1 /(M,) 

27 

1 +^i/(Mi) 


(109) 


2t 


where 


/(M,)=^2- 1 + /i +^2 (M, 2 — i)+(Mil_2) 2 

Expanding the right-hand side in powers of Mi 2 — 1 (assumed 
small), one obtains 

( 7 + 1 ) tan 0 = (M, 2 ~/ ) 3 / 2 [l-^^(M, 2 -l)+. . .] ( 110 ) 

Reverting this series and letting 0 «tan 0 and M l =M a>Sf 
one finds 


€ 


M^-l 


°° s [(7+l)0 ] 2/3 

For 7 = 1 . 4 , 


■ 2 1 /3 


1- 


I + IO7 


12 


(7T 


0 2/3 +O(0 4/3 ) 


(in) 


1+107 


12 


(+T 


1.176 


APPENDIX E 


VARIATION OF REDUCED MACH NUMBER BEHIND A PRANDTL-MEYER EXPANSION FROM M= 1 THROUGH AN 

ANGLE WITH 6 


The exact relation here is 


_ /7+I 


\ 7 — 1 


tan - 


. h - 1 

Vy+i v 


M pm 2 — 1 — tan _1 yiliV 3/ 2 — 1 (112) 


Expanding the right-hand side in terms of VM />Af 2 — 1, using 


_ 1 1 ^ ^ 

x = z_ "3 + 5'“T*’ ‘ 


one obtains 


( 7 + i)«=|(M fJ , 2 - ir s 1)n ( ii3) 

6 71 = 1 14- — 


where 




Reverting this series, the first few terms are 


Therefore, 


47 r3(7+D ft T /3 

"5W+i)l~^r- 01 


,J 3 + 0(04/3) | 

(114) 


i,pM - 


MV- 1 

[(7+1)0] 


2/3 


3\ 2/3 


,, 47 f3(7 + l)7 / 

+ 5( 7 +l)L 2 J 


0 2 «+o(0V3) (115) 


For 7 = 1 .4, 


47T3(7+l)-f 3 =1 Q97 
5(7+ 1)L 2 J 1Uy/ 
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